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ABSTRACT 



to the downstream direction of interstellar neutral 



The ESA meteoroid model predicts impacts of metc- 



Helium on i ts way through the heliosphere (Witte 



oroids in the mass range between 10~^® to 10*^ g on 
spacecraft surfaces. It covers heliocentric distances 
from 0.3 to 20 AU. Measurements of the dust detec- 
tor on board the highly successful joint ESA/NASA 
mission Ulysses have shown, that the flux of mete- 
oroids with masses between 10~^^ and 10~^^ g is, at 
least in the outer Solar System, dominated by in- 
terstellar dust grains that traverse the Solar System 
as it travels through the local interstellar cloud. We 
present a simple semi-analytic interstellar dust model 
that can easily be included in the ESA meteoroid 
model, together with a more precise determination 
of the flux direction of the interstellar dust stream. 
The model is based on the assumption that inter- 
stellar dust dynamics have two effects: solar gravi- 
tation and radiation pressure determines the spatial 
distribution, and Lorentz-interaction of the charged 
particles creates a temporal variation. 



1. INTRODUCTION 

As meteoroid impacts are a concern for highly sensi- 
tive instruments on board Earth orbiting satellites as 
well as interplanetary exploration probes, ESA has 
decided to develop a standard tool for the prediction 
of impact fluxes. In its first release version the tool 



(Staubach et al. , 1997) has been used to predict me- 
teoroid impact rates in order to assess its contribu- 
tion relative to the small debris population in Earth 
orbit. Due to the increase of data available from 
the in situ measurements of Ulysses and Galileo, as 
well as from the radar meteor measurements by the 
AMOR facility in New Zealand, ESA has decided to 
upgrade the tool. In the course of this update, an ex- 
tension of the interstellar dust module was proposed. 

The former implementation of the ESA meteoroid 
model includes an interstellar dust component that 
is modelled as a constant, mono-directional stream of 
particles moving towards a heliocentric ecliptic direc- 
tion of /3ecl — —5°, Aecl — 79°, with a flux density 
(integrated over all grain sizes) of 1.5 x 10~^m"^s~^ 
and a velocity of 26 km s~^. The direction is identical 



et al. , 1993 ). From the analysis of the Ulysses dust 
data it was found, that the downstream direction 
of interstellar dust is st atistically indistingu ishable 
from the gas direction ( Frisch et al. , 1999 ). The 
current implementation of the interstellar dust mod- 
ule assumes that the particles move on straight lines 
through the Solar System. This is a good approxi- 
mation for particles with masses around 3 x 10~^^kg, 
for which the forces of gravity and radiation pressure 
nearly cancel each other. The ratio (3 of radiation 
pressure to gravity is equal unity for these particles 
(see figure p. 
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Figure 1. Ratio (3 of radiation pressure force to 
gravity on dust grains in the Solar System as a func- 
tion of grain mass m. The dashed line indicates 
f3 ~ I, for which radiat ion pressure and gravity ex- 
actly cancel each other ( Gustafson, 1994 )- 



2. SPATIAL DISTRIBUTION OF 
INTERSTELLAR DUST IN THE 
SOLAR SYSTEM 

For values of f3 other than unity, the particles move 
on Kepler hyperbola (figure ||). As a consequence, 
the apparent stream direction and local concentra- 
tion depends on the location of the detector in the 



solar system (Landgraf et al. , 1999) 
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Figure 2. Spatial distribution and motion of inter- 
stellar dust grains with (3^1 m a 40 x 40 AU region 
around the Sun in the ecliptic plane. The gray-scale 
contours indicate the spatial concentration and the 
solid lines show the hyperbolic paths of grains with 
selected initial conditions. Diagram (a) shows the 
situation for grains with (3 = 0.5, and diagram (b) 
for grains with (3 = 1.5. 



TEMPORAL MODULATION OF THE 
LOCAL INTERSTELLAR FLUX 



a Lorentz force acts on the grains. The direction this 
force depends on the phase of the solar cycle. Dur- 
ing the period from 1991 to 2002 it is was directed 
such that interstellar grains are deflected away from 
the solar equator. Consequently the local interstel- 
lar dust concentration is significantly reduced. While 
proof of the electromagnetic hypothesis still pe nding , 
we employ the model described by Landgraf ( ^OOC ) 
to calculate the temporal flux variation for various 
grain sizes (and thus charge-to- mass ratios). Fig- 
ure 1^ shows the normalised spatial concentration of 
grains of various sizes as a function of time over one 
full solar cycle from 1991 to 2013. For flux predic- 
tions after 2013 we assume a perfect periodicity of 
the solar cycle. 
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Figure 3. Temporal variation of the flux of inter- 
stellar dust in the solar system (inside lOAUj over a 
full 22 year solar cycle for various grain masses. As 
a reference, grains with (3=1 and charge Q = C 
are given (thick, solid line) that are effectively not 
affected by gravity, radiation pressure, and Lorentz- 
force. 



Because of the differential susceptibility of small and 
big grains to the electromagnetic deflection described 
above, the grain mass distribution far from the Sun 
can not me measured directly inside the solar system. 
We assume that the grain masses are distributed ac- 
cording to 



Ulysses measured a variation of a factor of 3 in the 
flux of interstellar dust grains in the Solar System 
(Landgraf et al. , 2000). This was attributed to 



the electromagnetic interaction of the charged grains 
with the solar wind m agnetic flcld tha t varies with 
the 22 year solar cycle ( Landgraf, 2000] ). Interstellar 
dust grains are charged mainly due to the photo- 
effect caused by solar UV photons. Because of the 
magnetic field that is frozen in the radially expand- 
ing solar wind (speed between 300 and 700 km s^^). 



n{m)dm = nQUi ^^'^^dm. 



(1) 



We implement the temporal variation of the dust 
concentration in the interstellar dust module of the 
ESA meteoroid model by modulating the grain mass 
distribution over time. 



4. INTERSTELLAR DUST STREAM 
DIRECTION 

The interstellar dust stream direction has been deter- 
mined from Ulysses d ata collected within four weeks 
of the Jupiter fly-by ( |Frisch et al. , 19991) . With the 
now much extended database, a more accurate de- 
termination is possible. 
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Figure 4- Best-fit analysis of the Ulysses data 
between Jupiter fly-by and south polar pass. The con- 
tour levels (from light to dark) indicate 10%, 5%, and 
1% confidence intervals. Crosses mark local minima 
of the function. The lower right cross represents 
the global minimum (99% confidence level), the sec- 
ond local minimum represents a fit with a 98% con- 
fidence level. 

Figure ^ shows the fit-levels of an interstellar dust 
directional model compared with data provided by 
the Ulysses dust experiment in the time range from 
the Jupiter fly-by in February 1992 to January 1994, 
just before the first pass of Ulysses over the south 
pole of the Sun. For the analysis we use only im- 
pacts with a mass (also measured by the detector) 
between 10"^** kg and 10~^^ kg. For these grains 
the solar radiation pressure and the gravity of the 
Sun are equal (figure 0). Thus, it can be assumed 
that such grains are not deviated from their orig- 
inal direction of motion. Moreover, only impacts, 
which are the most reliably detected, are used for 
the comparison with the model. The interstellar 
dust flow model can be characterised by the vector 
(Aecl,/?ecl,o,o-/3), where Aecl and /3ecl,o are re- 
spectively the ecliptic longitude and latitude of the 
dust flow direction, and ap the width of the direc- 
tional distribution. From this model the total num- 
ber of dust impacts for each spacecraft rotation an- 
gle interval is calculated. The expected number of 
impacts per rotation angle, derived from the model 
and the measured number of impacts per rotation 
angle interval have been compared by the sta- 
tistical method. Its minimum, presented in figure 
^, gives the model parameters providing the best fit 
to the data. The best-fit parameters for the main 



downstream direction of interstellar dust in the So- 
lar System is hehocentric longitude Aecl — 110° and 
hehocentric latitude /3ecl,o = —30°. Also a broad- 
ening of the directionality by ap — 20° was found to 
provide the best fit. 

It can be seen from figure ^ that the best-fit flow di- 
rection of interstellar dust is not determined unam- 
biguously. Another good fit is achieved for Aecl = 
30° and /3ecl,o = —10°. Analysis of the Ulysses data 
is ongoing in order to resolve this ambiguity. 
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